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nonrelativist ic  method  has  also  been  found. 


Is  C  an  adiabatic  invariant  for 


relativistic  particles  slowing  down?  Yes.  A  modification  of  the  method  of  Alfven 

I  4 

and  Falthammar  that  shows  this  is  presented  in  the  next  section.  This  is  fol¬ 
lowed  by  a  rigorous,  longer  method  that  r  esults  in  the  time  derivative  of  O.  Other 
adiabatic  invariants  are  then  examined  briefly.  A  condensed  version  of  this  report 
has  been  submitted  to  a  journal.  ^ 


2.  SIMPLE  DERIVATION 


At  any  time  t  choose  the  origin  of  a  cylindrical  coordinate  system  at  the 

center  of  gyration  of  the  particle  of  charge  e  and  rest  mass  m  with  the  z  axis  in 

the  direction  of  B.  Assume  a  constant  magnetic  field  gradient  c — .  From 
BBx  9Br  R  9B  3z 

V-  B  =  0,  neglecting  — and  — -a,  BD  =  -  v  ^ .  v,  =  v.;  ev,  and  thus  ep  are 
*  6  9^  R  2  9z  i  <f>  <f>  ri 

negative.  Let 


*  _  dX 
X  '  dt 


P  =  Plf  +  Pi 


2,2 

1  -  v  /c 


r  J  2  4  2  2  2  J  2  4  2  2  2  2 

E  -  ^  m  c  +  p  c  =  myc  =  y  me  +  p|(c  +  p^c 


E  =  EP£_  =  m73vv  *  v(|  p„  +  vipi 


•  3  .  my  c  3« 

p  =  my  vv  — *  x  =  my  v. 


The  component  of  force  in  the  z  direction  is 


vn  eviR  a-R  v  |i 

P||  =  ■ev1  Br  +  p  —  =  —3—  ^  +  p  — 


V1  Pi  3B  ,  3  vll  • 

”  TT  B  dsT  +  my  —v  <8) 


5.  Dubs,  C.W.  (1980)  Adiabatic  invariants  for  a  charged  particle  slowing  down, 
submitted  to  Journal  of  Geophysical  Research. 


since  R  nearly  -  -  — g ,  the  radius  of  gyration.  Substituting  this  into  the  expression 
fur  6: 


_  3  •  Vipi  dz  913  ,  3  vll  •  ,  i 

vv  : ar  —  v  +  vipi- 

2  2 

3  v  '  vn  .  .  v  p  . 

— —  v  +  vipi  '  TR-®  0  • 


(9) 


(10) 


Multiplying  by 


ViPl 


2-1-V  +  —  p-i-B-0  . 
v  p  1  1  B 


(ID 


Substituting  p^  m^v  sin  a: 


*> 


iZ 

.  •> 
v  +  — 
>v 

d(->v)  4  2 

d  sin  a 

1 

dB 

2  d  sin  o 

1 

dB 

0 

V 

dt  sin  a 

dt 

B 

dt 

sin  o  dt 

B 

dt 

(12) 

1 

d_ 

dt 

In 

9 

sinw  o 

B  J 

0  , 

(13) 

sin' 

) 

“  a 

const 

ant  . 

(14) 

B 

X  KKiOHOl  S  PKRIV  ATIOIN 


3.1  Simple  Mapm'lie  Field 

Choose  a  cylindrical  coordinate  system  as  above.  Assume  a  vector  potential 


(15) 


bet  bB  be  the  magnetic  field.  The  normalized  magnetic  field  then  is 


B 


k  r 


(in) 


(II  is  normalized  from  here  through  Kq.  ((13). )  The  constants  b,  h,  and  H  are  seen 
to  be  the  field  at  t lit'  origin  and  the  B^  doubling  distances  respectively. 
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L  =  -  m  c  ^  1  -  v^/c^  +  eA  '  v  ,  v=RR+$R$+zz 


d  / 3L  \  9L  ~ 

*  UT)  '  ^  Q‘  ' 


(17) 

(18) 


represents  the  generalized  component  of  the  drag  force  f(v)v.  (See  Section  3.2.) 


-jp  =  m^vv  =  f  »  v  =  f( v) v  • 


(19) 


So, 


f(v)  =  my  °v 


(20) 


(The  latter  equals  p,  as  it  must. )  Only  the  R  and  z  Lagrange  equations  are 
needed: 


^  ImyR] 


Tt  l®**! 


ebR^J  , 

2h 


3.  z 
my  v  — 
v 


(21a) 


(21b) 


which  yield 

(22a) 

(22b) 

The  equations  are  now  partially  normalized  with  e  =  -j—  ,  E  =  ,  R  =  R^p, 

z  =  Rq£,  v  =  Rqu,  where  Rq  is  the  initial  value  of  R  and  of  the  gyroradius.  (Note 
that  E  from  here  to  the  end  of  Section  3.  1  means  this  instead  of  energy.  ) 


R  =  R  4>  + 


e  b  R  j>  .  v  +  vR 

my  '  '  v 


e  b  R  4>  v_z 

2myh  v 


(  )  * 

p  =  p 


^  1  +  jj-  +  Hy  =  (1  +  €  S  +  Ep) 

;2  +  ( )  +  if 

y  my  u 


(23) 


(24a) 
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1 


€ e  b p  “  o  4  it 
*Tm^  u 


l)  -  h 


dC  ci_ 

lit  dt 


r  i 

- 

- 

1  1  -  (’|1S“  o 

.1 

l 

(>.*•  ii>:! 

L  11 . J 

<it 

j 

B 

uL>  B3 

(24  b) 


(2  0 


l)  -  A  ^  *  z  (  >  . 


Ci'O 


K“  (  )-  4  , 


(27' 


u  K  p  4  o  p  6  *■  y.S 


(28) 


u  •  B  ?  (  )  -  ^  p  p  . 


(2P) 


If.,  2  2  1 

-1/2 

r  1 

9  9  ' 

r;1/2l 

1 

?()--§■  p  p 

u 

<r  +  V1 

1  J 

B 


)  (c  S  +  K  p )  +  p  p 


2  3 

u  B 


n ) -f  pp 


f  (  >  +  £  (e  S  +  Kp)  -  (p“  +  pp) 


2ti 

3  t,3 

u  B 


“ 

2  3 

■ 

2 

r  .  2  i 

?< )  -  |pp 

+  2„5 

u  B 

co  -  f  PP 

(  )  (e  ?  +  Ep)  +  ypp 

(30) 


2  3 

Lot  M  be  u~  B  times  the  middle  two  of  these  four  terms.  Use  the  two  equations  of 
motion  to  eliminate  p  and  S. 


M  =  -  2 


?  (  ) 


P  P 


L  2m-v 


(  ) 


eebp^  ,  ,  eppit 
2m>  2u 


Lk  M  +  £  P  k  H 
u  w  2u 


(31a) 


c(u2  -  3  ?2)  -  2Ep  ? 


(3  lb) 
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\  'r  that  all  of  !li»*  relat  i vist  ie  terms  and  all  of  the  terms  containing  the  •  hange  in 
ti.e  particle's  speed  cancel.  So,  D  is  the  same  with  or  without  slowing  down  and 
with  or  without  being  derived  relat ivist ically.  The  result  then  is: 


D 


_L_  J.Zsl 

u2B5  4 


p(u~  -f  £“)(!  +  e  £  +  Ep)“  +  e  p  r^(  1  +  e  £  +  Ep) 


€  w  2  .  9-9  *,  9 

+  “x*p  p(pw y  -  o 


+  Kp 


-  (p  ~  +  p  ~  0 w)( 1  +  e  £  +  E  p  )2 


-2ep^(l+e(+Ep)2  + jp2(2p2 


2  •,  2 


p^"  -  3  r2xi  +  £  r  +  Ep) 


(32) 


2  3  4 

Note  that,  for  E  =  0  (H  =  oo),  D  contains  terms  only  in  e  ,  e  and  e  divided 
by  B  ,  and,  for  e  -  0  (h  -  oo),  D  contains  terms  only  in  E,  E2,  E2,  and  E^ 
divided  by  u“  Alternatively,  D  -  0  (e“  E°)  +  0  (c°E). 


3.2  General  Magnetic  Field 

Which  of  these  results  hold  for  any  magnetic  field?  Repeat  the  last  section, 
but  with 


A  -  Hp  0  +  e  ba* (H,  0,  z) 


a  r  a'  +  VT 


(33) 


where  a1  is  an  arbitrary  function  of  R,  0,  and  z  except  that  vXa'  vanishes  at 
R  ~  z  =  0.  Then  bB  V  X  A  is  completely  arbitrary,  regardless  of  what  function  is 
chosen  for  the  scalar  V.  To  simplify  the  algebra,  choose 

v  =  -  /  a}?  dR  ,  (34) 

holding  0  and  z.  constant.  Then  =  0,  and  the  normalized  field  is 

3(Ra , ) 


K  ^rxA-z  +  €VXaT=z  +  e 


3a  0a . 


~  f  1  Wcl  \  I  on 

r(b  ^-ai4)  +$ 


R  00 


da 


OR 


-  1 


R  OR 


h  R  e  +  0  e  jT  +  z( !  +  c  d?) 


(3  5) 
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nnHmn  V  ft  +  ‘in*,  — 


R  ~  R  0 


2  *-  —  In  o  *  e(Ro  .:  -  '/ .:  )l 


/  0  V 


In  general,  the  equation  for  H  0  if*  needed. 


d_  /  9L  \  9h  +  Q 

dl  Vai  j  »•  * 


|  mTR2«  +  eb  +  eRaJ  =  R  i  ^  *)  +  m‘’3  '~T^ 


.  2.  |  d(Ra6> 

m7  R2  6  +  2nr>  R  R  i  +  mV  ^R“i  +  cb  RH+f  jt — 2~ 


/  9a,  da  \  <; 

eob(  aT  R*  +  W  z)  +  mV 


(  ta.  i) a  \  d(  Ra  , ) 

rt  •  7'1  O 


■2  R  6  +  -  ft  +  ^  t  TTa  h  0  * 


)6  /  dt 


Rd  -•  -2ft*  +  ££  l-R  +  f  0.)R  -  R-yi  +  ^ 


<L  *L]  SL+Q  . 

dtU*J  3z  z 


9a,  9  a 


^  (rn*>  z  +  obf  a^)  *bc  (  Ro  +  z  )  +  m*>  v 


■S  +  nn3  4  /.  +  ceb/  f  eb  ^  R  i  +  *T  *  I  +  "n  ^ 


9a ,  9a 


eb/f^  _i  +  !fz.  !!za  ^2 

'  e  my  l  9z  R<i>  dz  Z  '  8R  R  '  8$ 


z  =e  ^<-r^r  +  «V  +  TT  • 


^3Z  "I  V  z 
3z  ZJ  v 


<5G) 


(57) 


B2  =  e2(/32  +  /32)  +  (1  +  e(3z)2  . 


v  •  B  =  z  +  c(R/3r  +  +  z(3z)  .  (58) 


D  =  {(B2)'1/2  -  [z  +  e(R0R  +  R*fy+  z0z)2  v'2  (B2)'3/2  j 

(t  ^egz)egz  +  e2(3R  3R  +  ^  _  2[  1  cCft  3R  +  R R +  z  jg 


v2  B3 


2[  JfSll  +€0z>*tZ(iR*cRiB4}  ,  j  2(_2  y) 

-  - - * —  -  *■ 


2  -p3 

v  B 


TZF 

v  B 


[  )"  (-3)1(1  +<3z>e&z  +  e2(0R3R  +  0*/y} 

JZ3 - 

v  B 


(59) 


where  [  ]  is  z  +  e(R  0R  +  R  4>  0^  +  z/g.  Using  Eqs.  (47),  (52),  and  (57),  -  v2  B3 
times  the  third  and  fourth  terms  is 


-t  i 


(c  5§-<-R*0r  + A04XI  +e3z)+  ^<l  +  e<y +  eR*20f 


+  e  ~  <R$0  R  +  eRi0  R0  _  -  e  z0  R3  J  +  e  ^  3R  -  2  c  R  i  0 


vR 


+  C^<-6V”%Ve*Vz)+e 


V  -«  (^R^V^z)} 


vz  /vR 


=  -[  ]  e(Ri2|3R  -  2Ri3,}  . 


(60) 
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r 


Note  two  things:  First,  the  v  terms  cancel,  so  I)  is  unchanged  by  I'nnrs  hi  the 
tv*  direction  and  thus  of  any  processes  producing  only  slowing  down.  Second,  the 
terms  containing  •>  cancel,  so  D  is  the  same  with  a  relativistic  derivation  as  it  is 
with  a  nonrelativistir  one.  So 

2/ 


u  -  -  •  1  I— 


B 


2  [z  f  c(R  dR  +  R  i  +  z  t7\e  {R  +  H  cT  dR  -  R  4  *■  R  ©  *  z  ( 


v2B3 


3{y  +  H  0  J  f  )“{<!  f  f  j  )f,i  *■  c“(dp  dl?  *•  o’  S  >  } 

_ K _  <j)  v.  ^  /  Z _ 1  v  K  0  0' 


Hdl{  r  1UMh  -  %} 

I)  -  2c  /.  - - - ^ - £ - 2  f  o  <c~)  . 


Dig) 


Thus,  in  general,  1)  is  of  first  order  in  c. 

Perhaps  the  simplest  way  to  see  that  I)  is  independent  of  slowing  down  process¬ 
es  is  as  follows. 


f) 


B 


B 

1  2 


cjv 

dt 


(<;:0 


Since  IT  (unnormalized  from  in*re  on)  is  a  function  only  of  K,  <j>,  and  the  only 

dv 

thing  in  this  expression  which  could  depend  on  slowing  down  is  -rr  . 


<Jl } 

dt 


HI-}  V 


dv 

dt 


d(i  v) 
dt 


e  v  X  B  +  f  v 


CM) 


Dotting  v  into  this: 


r 


m 


dfrv) 

dt 


<•>;>) 
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civ  * 

since  yjy*  must  be  perpendicular  to  v. 


So 


dv  e  v*  X  B 
dt  myv 


<m>) 


the  same  as  it  would  without  slowing  down.  Since  this  is  the  only  term  in  I)  con¬ 
taining  7  but  it  contributes  nothing  to  D,  the  latter  is  seen  again  to  be  independent 
of  whether  or  not  it  is  derived  relativistically. 

If  there  be  a  component  of  drag  force  perpendicular  to  v,  write  the  drag  fort  e 
as  f  =  fyv  +  f^n,  n  *  v  =  0.  Carrying  out  the  same  procedure, 


dv 

dt 


e  v  X  B  +  f  n 
_ n 

myv 


(67) 


and 


B  • 


dv 

dt 


f 

n 


myv 


B  •  n 


instead  of  zero  in  the  expression  for  D. 


(68) 


4.  OTHER  ADIABATIC  INVARIANTS 

More  generally,  consider  a  purely  field-geometric  quantity,  g(B,p),  p  =  cos  a. 


By  the  same  reasoning  (p.  14,  last  paragraph),  g  is  seen  to  be  independent  of 
slowing  down  and  of  relativity.  Further,  consider 

s' 

m 

G  =  f  g(B,  jj)  ds  . 
s 

m 

where  s  is  the  distance  along  a  field  line  quite  close  to  the  locus  of  guiding  centers, 

and  sm  and  are  mirror  point  values.  For  C  =  -g—  -  const.,  Bm  the  mirror 

m 

point  value  of  B,  s  and  s'  are  constant.  Then 
m  m 
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1  ♦ 


A  treatment  like  that  in  Section  *.  1  leads  tt 


e  b  H  “  $  / 

?  - — y —  -  ft 
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Suppose  that  at  t  0:  v  v  ,  z  0,  /  /  ,  and  H  H  ,  the  radius  of  gyration. 
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Integrat  ing: 


mvNC 


M 


Jin  H 


ny 


2  7 Jkv  y/Z7 


All  three  of  these  i|uanf  iftes  are  seen  f.<  var\  w  ith  rime.  No,  even  though  thrv  .irc 
adiabatic  invariants  with  parallel  n  mpi  ‘tiers!  >  .  .f  fi.n  (<  the\  ;it  r  not  with  slow  mil' 
down.  K  wil  be  essentially  invariant  f.  it**  \  but  an  t  «  ases  oxfiotirnt  la  I  |v  w  ill 
t  at  large  times  unless  k  i  is  a  haP’  mt«  •  ‘lines 


X  CONCH  SIONS 

rhe  first  adiabatic  invariant,  :>  J*  Ik  i  n  >*  inv.i'  mill  if  'In-  particle  is  .-l.-winr 
down.  The  eloselv  related  <  ftintif  1 1  \  1  —  L'  .  lo\\<  \*  it  is  an  adl.ibaMi  imariant 

f«>r  anv  energy  particle  in  anv  magnetic  field,  r\«  n  ,\ith  an  additional  force  -d‘  tin 
form  f(v)v.  This  is  not  surprising  ^im«  a  f-s.e  m  the  t  v*  di  n-i  t  n  >n  should  n<  >t 
change  i,  and  thus  nof  (  ,  a  1 1  h<  *ui»h  d  w  .  add  •  a  use  the  *  rtit  t-r  of  g  \  i  at  i<  *n  to  i *  c  >v* 
across  the  field  (toward  the  particle  if  it  is  brine  slowed  d«*wnl  since  the  rviai  ndm.s 
is  proportional  to  v  .  (  is  not  invar  iant  with  a  parallel  «  ot  ipononl  of  rln  ti  n  field 
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t 


ff 
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present.  D  h  b-^p-,  b  a  constu*  ,  i?  shown  to  be  independent  , ,f  f  and  thus  of  slow  - 

ing  down  for  all  magnetic  fields.  It  m  also  shown  to  be  the  same  whether  derived 

rel  a .  i  v  i  -0  it;i  11  y  a  not.  Ip  y.  iernl,  D  is  /<*.-,  ,  ip.*  first  .  r*  f  !*  ♦  /radjerv 

For  the  simplest  *  magnetic  field  with  a  paial:  1  gradient  of  B^,  t ,  and  perpon- 

dicular  gradient  of  B  ,  —  ,  D  0  (c^E°)  +  0  (e°E),  where  r  H  'h  E  R  H 
z  11  o'  o' 

R  is  the  initial  gyroradius,  and  B  bz  at  the  origin. 

The  adiabatic  invariance  of  C  with  a  force  on  the  particle  in  the  tv  direction, 
among  other  cases,  is  applicable  to  trapped  and  auroral  protons  between  nonfor¬ 
ward  scatterings.  Examples  of  processes  which  may  cause  charged  particles  to 
slow  down  without  large  angle  scattering  are  elastic  and  inelastic  scattering,  ioni¬ 
zation,  Bremsstrahlung,  and  Cerenkov  radiation. 

In  general,  field  geometric  quantities:  C,  half  bounce  path  length,  I,  and  $ 
remain  adiabaticallv  invariant  with  slowing  down,  but  M,  J,  and  K  do  not. 
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